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We use the light-front Hamiltonian of transverse lattice 
gauge theory to compute from first principles the glueball 
spectrum and light-front wavefunctions in the leading order 
of the 1/N C colour expansion. We find ++ , 2 ++ , and l + ~ 
glueballs having masses consistent with N c = 3 data available 
from Euclidean lattice path integral methods. The wavefunc- 
tions exhibit a light-front constituent gluon structure. 

PACS numbers: 11.15.Tk, ll.10.Jj, 12.38.Bx, 13.85.Ni 



I. INTRODUCTION 

There is growing experimental evidence that glueballs, 
boundstates of gluons in the SU(3) gauge theory Quan- 
tum Chromodynamics (QCD), have been discovered in 
the mass range 1.5-1.7 GcV [|l| | . But the confinement 
feature of QCD — interactions grow stronger as the en- 
ergy of a process decreases — complicates any first prin- 
ciples calculation of the boundstate problem. To date, 
the most successful boundstate calculations use a Eu- 
clidean spacetime lattice and simulate the path integral 
by Monte Carlo methods. These difficult calculations are 
now roughly consistent with the experimental signatures 
of glueballs although much detail remains to be 

clarified. Therefore, it is important to have some inde- 
pendent method of calculating the properties of glueballs 
from first principles in gauge theory. In this letter, we 
present our first results on this problem using an effective 
light-front Hamiltonian quantisation (canonical quanti- 
sation on a null- plane in spacetime). This is the trans- 
verse lattice method, suggested originally by Bardeen and 
Pearson (5|(| , which we have developed to the extent that 
quantitative calculations are now feasible [[?]||] . 

Although this work is ostensibly about glueballs, we 
have a more general motivation for developing the light- 
front Hamiltonian formulation of gauge theory. A light- 
front Hamiltonian has Lorentz-frame-independent wave- 
functions. Together with the simplicity of the light-front 
vacuum, this leads to a field-theoretic realisation of the 
parton model for hadrons on which so much understand- 
ing is based. The rich phenomenology in hadronic and 
nuclear physics that would follow from knowledge of the 
light-front wavefunctions is surveyed in the lectures of 
Brodsky ||. We will use the glueball problem in QCD, 
which is especially difficult computationally, as a quanti- 
tative test of our light-front formalism. 



A detailed account of our methods and various quan- 
titative tests that we have performed, mostly in 2 + 1 di- 
mensions, can be found in Refs. but we will briefly 
review the salient points below. We will work in the 
leading order of the 1 /^-expansion of SU(N C ) gauge 
theory, which omits the l/N 2 -suppressed glue configura- 
tions. It also removes "sea" quarks from our theory, an 
approximation also used in the Euclidean lattice calcu- 
lations. We work on a coarse transverse lattice, using 
an effective potential tuned to minimize discretisation 
errors. The groundstate J vc = ++ glueball mass is 
M = 3.3 ± 0.2^, where a « 0.1936GeV 2 is the string 
tension, and is consistent with SU(2) and SU(3) results 
available from the Euclidean lattice. Components of the 
2 ++ and l" 1 glueballs which are behaving covariantly 
are also consistent. The glueball wavefunctions that we 
obtain are new. 



II. TRANSVERSE LATTICE FORMULATION 

In 3 + 1 spacetime dimensions we use a square lattice 
of spacing a in the 'transverse' directions x = {x 1 ^ 2 }, 
and a continuum in the {x°,x } directions. Our action 
is 

S= J dx°dx 3 ^2 yTi^D a M r (x) ^D a M r (x)) t | 

where r, s € {1, 2}; a, S {0, 3}; and 

D a M r (x) = (d a + iA a (x)) M r (x) - iAf r (x)A Q (x + of) . 

(2.2) 

A a (x) is the gauge potential at transverse site x, while 
M,.(x) £ GL(N C , C) is a link variable from x to x + ar 
which carries colour flux between sites. Our use of non- 
compact variables at finite a is the basis of the colour- 
dielectr ic fo rmulation of lattice gauge theory |10|. The 
action (|2.l[ ) is gauge invariant at each transverse lattice 
site under 



A Q (x)->[/(x)A Q (x)[/t(x) 
M r (x) -> C/(x)M r (x)C/ t (x 



H(<9 Q t/(x))[/t( x ) (2.3) 
ar) (2.4) 



for U E SU(N C ). The dimensionful coupling G 2 (a) is 
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a — > 0, where g 2 is the usual continuum gauge coupling. 
We include a gauge-invariant effective potential V^[M] 
consisting of Wilson flux loops on the transverse lattice. 
The exact form we use is given below. Part of its job 
is to ensure that M is forced to be a unitary matrix as 
a — ► 0, so that continuum QCD is recovered a la Wilson 
Jj~i] . More generally, we choose it to minimize violations 
of Lorentz covariance at finite a. The idea is that we 
tune the couplings of V along an approximation to the 
Lorentz covariant scaling trajectory which leads to con- 
tinuum QCD as a — > 0. Then for low energy observables, 
such as boundstate masses, one can work at finite a. 

We perform light-front quantisation, treating x + = 
(x°+x 3 )/ y/2 as time, and pick the light-front gauge A_ = 
(A - A 3 )/V2 = 0. In this gauge A+ = (A + A 3 )/V2 
is a constrained field which we eliminate by solving its 
equation of motion classically. Working in the leading 
1/N C approximation has some simplifying consequences. 
In the transverse rest frame P = 0, Eguchi-Kawai reduc- 
tion occurs at large a, so that the argument x may be 
dropped JtJ. This results in a light-front Hamiltonian of 
the form 



r , dx --^ T J J+1 Lj + 
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4/V ' Tr|MiM^| Tr|M 2 Af|} , 



a 2 N 2 

J+ = i(M r d- Ml + Ml d- M r ) . 



(2.5) 
(2.6) 



This Hamiltonian includes all Eguchi-Kawai-reduced 
Wilson loops on the transverse lattice up to fourth order 
in M; this is our approximation to the effective poten- 
tial V. Generalisation to higher orders is straightforward. 
The Hamiltonian is diagonalised in a parton Fock space of 
M r , at fixed total momentum P + , built from its Fourier 
modes at x + = 



M r 



1 



dk+ 



(a- r 



e~ lk+x +aU ik+x 



(2.7) 



However, only closed flux loops, invariant under residual 
^^-independent gauge transformations U , have finite en- 
ergy. These loops propagate without splitting or joining 
at leading 1/N C order. 



Even with a finite lattice spacing a, the Hamiltonian 
THI) is not yet completely regulated. We impose periodic 
boundary conditions in x~ and a cut-off on the maximum 
number of link partons a r in the Fock space. For each 
choice of couplings in the Hamiltonian, we remove these 
latter cut-offs by extrapolation, following a continuum 
improvement procedure. 



III. TUNING THE EFFECTIVE POTENTIAL 

To test for Lorentz covariance, we measure the dis- 
persion relations of glueball boundstates as a function 
of couplings in V. This requires us to work at P ^ 0, 
which can be achieved by adding appropriate phase fac- 
tors to matrix elements of P~ above ||. Using G 2 N C to 
set the dimensionful scale, the dispersion relation of each 
glueball can be written 



2P+P- = G 2 N C (Ml + Mj a 2 \P\ 



+2M 2 1 a 2 P 1 P 2 + 0(a i \P\ 4 [ 



(3.1) 



For each glueball in the light spectrum and for fixed fi 2 , 
we tune the couplings A 2 ; and j3 so that, as far as possible, 



a 2 G 2 N c M\ = c 2 on = l 
a 2 G 2 N c {M\+M\) = c 2 oS = l 



(3.2) 
(3.3) 



Con is the speed of light in direction x = (1,0); c s is the 
speed of light in the direction x = (1, 1). In practice, we 
apply a % 2 -test over a range of light glueballs to optimize 
isotropy of the speed of light. In this way we search for 
a Lorentz covariant scaling trajectory parameterized by 
/i 2 , which may later be related to a. 

It is necessary to know the dimensionless combination 
a 2 G 2 N c for this procedure. This can be deduced from 
two measurements of the string tension a. The first mea- 
sures the mass squared of winding strings in the trans- 
verse lattice directions x, fit to n 2 a 2 cr| n for winding num- 
ber n. The second measures the heavy quark potential 
in the continuous x 3 direction fit to o^i? for sep- 
aration R. Demanding <jt = &l = & and expressing 
eigenvalues in units of G 2 N C , we may eliminate a and 
deduce a 2 G 2 N c . We may also deduce G 2 N C , and hence 
all dimensionful quantities, in units of a. 

To further help optimize lattice discretisation errors, 
we include in the % 2 test a measure of deviations from ro- 
tational invariance of the heavy-quark potential, includ- 
ing deviations from <jt = cl and deviations from Lorentz 
multiplet structure in the glueball spectrum. <jl/G 2 N c 
is the only basic quantity which is poorly determined by 
our current approximations, and our results were partic- 
ularly sensitive to the variance assignment in the x 2 test 
for this datum (though still compatible with the system- 
atic error estimate we make below). 

According to the colour-dielectric picture, large [i will 
correspond to large a where discretisation errors grow, 
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while small /i corresponds to intermediate spacing a. Suf- 
ficiently small a produces a non-trivial light-front vac- 
uum, signalled by the appearance of tachyons; it would be 
difficult to construct an effective potential in this regime. 
Instead, we use an effective potential at intermediate a 
on a trivial light-front vacuum. In 2 + f dimensions the 
above procedure accurately identified a Lorentz covariant 
scaling trajectory at intermediate a The results we 
present below are based on a similar trajectory which we 
have found in 3 + f dimensions, though further work is 
needed to reduce systematic errors. 



IV. RESULTS 

In Table | we show results for various components 
of the lightest glueballs at a transverse lattice spac- 
ing a = \.\j\J~o « 0.49 fm, at the overall x 2 min- 
imum for this a. In order to establish the scaling 
trajectory, we have performed calculations at other 
lattice spacings: At smaller lattice spacings tachyons 
appear, indicating breakdown of the trivial vacuum 
regime. For larger lattice spacings the couplings evolve 
gradually, and the minimum \ 2 value grows. For 
the groundstate ++ we find approximate scaling and 
maintenance of Lorentz covariance (a/i/tr, M/^/a, c) = 
{(LI, 3.3, 1.13), (1.2, 3.1, 1.18), (1.3, 3.5, 1.19)}. The er- 
rors shown on glueball masses in Table | are estimates 
of the error in removing the periodic boundary condi- 
tion on x~ and the cut-off on the maximum number link 
variables M r in a wavefunction; these errors are under 
control. Only for the ++ mass are we confident enough 
to estimate the finite-a systematic error. The Lorentz 
covariant scaling trajectory will need to be determined 
more accurately before similar statements can be made 
for the higher states. In 2 + 1 dimensions, we found 
that when the deduced speed of light becomes isotropic 
to within about 10%, the corresponding glueball mass is 
also accurate at this level; this also seems to be the case 
in 3+1 dimensions. 
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TABLE I. 

masses M 2 - 



The glueball multiplet components showing 
2P + P~(P = 0) and c from the dispersion re- 
lation (1.1). Vi is the symmetry under reflections x 1 — > — x 1 ; 



"Im" indicates c < and "— " indicates the quantity has not 
been measured. The J z = ±1 states are exactly degenerate. 
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FIG. 1. The variation of glueball masses with N c . ELMC 
predictions are continuum ones for N c = 2, 3, and fixed lattice 
spacings (higher points are j3 = 11.1 and lower are /3 = 10.9) 
for N c = 4 N c = oo light- front predictions are from 

components behaving covariantly in Table [l| including the 
total error estimate for the ++ . The dotted lines are to 
guide the eye, and correspond to the expectation of leading 
linear dependence on 1/A/?. 

We now compare our results with spectra from con- 
ventional Euclidean lattice Monte Carlo (ELMC) simu- 
lations. See Fig. |l]. (Recently, large N c glueball mass 
ratios have been calculated using a conjectured duality 
with supergravity fll3| ; we will not discuss these results 
since their errors are not quantified.) 
0++ 

Based on violations of Lorentz covariance and variation 
of M along the scaling trajectory, we include an estimate 
of systematic error for this state, M. = 3.3 ± Q.2y/a. A 
continuum extrapolation of ST/ (3) ELMC data has been 
made in Ref. Q, M = 3.65(15)^ (statistical error ) . 
Furthermore, measurements for SU(2) and S[/(4) |15| | 
indicate a slight downward drift with decreasing 
2++ 

Not all components of this state are behaving covariantly. 
The "Pi-odd combinations of J z — ±1 and J z = ±2 have 
large lattice artifacts and cannot be trusted. As a re- 
sult, we are unsure of the correct parity assignment for 
the J z = ±1 but assume here that it belongs to the 2 ++ 
multiplet. The 'Pi -even J z = and J z = ±2 combi- 
nation are nearly degenerate and isotropic, pointing to a 
tensor mass near M » 4.5y/a. The SU(3) ELMC extrap- 
olation of Ref. |0 gives 5.15(23), consistent with more 
recent coarse anisotropic ELMC simulations [^6|. Again, 
there is a slight drift with 1/N%. 
1+- 

Here the Pi-odd components fair better. The J z = 
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component of the lightest spin 1 is behaving covariantly, 
and points to M. w 5.9y/a. The J z = ±1 compo- 
nent is split away as a consequence of its anisotroj 
An SU(3) value M = 6.6(6)y/a is given in Ref. 
while the more accurate measurement of Ref. ]l6| implies 
M = 6.30(4)V^- 

o-+ 

Our lightest candidate for this state appears around 
7 Ay/a though we have not measured its covariance. 
SU(3) ELMC work suggests a lighter mass j|. Barring 
an unusual large-iV c extrapolation, it seems likely that 
this state (which posed problems for early ELMC studies 
also) is poorly described with our current approxima- 
tions. 

As an example of what may be deduced from wavefunc- 
tions, in Figure ^ we plot the distribution of longitudinal 
momentum fraction x — k + /P + in the ++ glueball at 
a = \.\jyfa. Gd{x) is the probability distribution for a 
link parton a r (k + ) and becomes the usual gluon distri- 
bution function in the continuum limit a — > 0. At fi- 
nite a, Gd may be interpreted as a light- front constituent 
gluon distribution that evolves with a. We find that our 
0++ light-front wavefunction consists mostly of four such 
constituents; it does not obey a gluonium picture (two 
constituent gluons). 
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FIG. 2. The ++ momentum distribution function. 



V. CONCLUSIONS 

We have measured glueball masses in QCD, using the 
l/N c expansion and different quantisation, field vari- 
ables, regulators, and gauge fixing from the traditional 
Euclidean lattice path integral simulations. The results 
from both approaches are consistent for QCD without 
quarks, indicating a ++ glueball near 1.6 GeV, and 
1/N% corrections appear to be small for the ++ and 
2 ++ states. (In this latter regard, it would be interesting 
to see better data at SCI (4).) We found that light-front 
wavefunctions exhibit a constituent gluon structure, al- 
though not of the form one might naively guess. In order 
to improve the accuracy of the Lorentz covariant scal- 
ing trajectory, future work includes choosing other ob- 
servables to set the scale and incorporating higher or- 



der effective interactions. Finally, we note that the re- 
introduction of quarks will be necessary for a detailed 
comparison with experimental glueball candidates, which 
mix with nearby scalar quarkonia 
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